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The differential intensities of Cosmic Rays at Earth were calculated using a 2D 
stochastic Montecarlo diffusion code and compared with observation data. We 
evaluated the effect of stretched and compressed heliospheres on the Cosmic 
Ray intensities at the Earth. This was studied introducing a dependence of the 
diffusion parameter on the heliospherical size. Then, we found that the opti- 
mum value of the heliospherical radius better accounting for experimental data. 
We also found that the obtained values depends on solar activity. Our results 
are compatible with Voyager observations and with models of heliospherical 
size modulation. 
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1. The 2D Model of the Heliosphere 

HelMod Code^ solves the bi-dimensional Parker's particle transport equa- 
tion^. A Monte Carlo technique is applied on a set of Stochastic Differential 
Equations (SDEs) fully equivalent to the Parker's equation'^. The model 
takes into account particle drift effects and latitudinal dependence of the 
solar wind speed and of the Interplanetary Magnetic Field (IMF). It is de- 
scribed in details in Ref. 1. In the model, the IMF from Parker"' is modified 
introducing a small latitudinal components as described in Ref. 5. For pe- 
riods of low solar activity, we take a solar wind speed gradually increasing 
from the Earth position up to a maximum value near the heliospherical 
poles (~ 760km/s)^. For periods approaching the solar maximum we as- 
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sume a solar wind speed independent on the latitude. 

The symmetric part of the diffusion tensor, in a reference frame with one 
axis aligned with the Parker's magnetic-field, is purely diagonal containing 
transverse {K±0 and K±r) and parallel (^||) components''. The diffusion 
coefficients are given by* 



K±r = Pk , 

K±0 = l{6) pk K\\ , 

where j3 = v/c,v the particle velocity and c the speed of light; the diffusion 

parameter accounts for the dependence on the solar activity; i?® is 
the measured value of IMF at the Earth position - typically « 5 nT, but 
changing with time - obtained from Ref. 9; B is the magnitude of the large 
scale IMF as a function of heliocentric coordinates; finally, the term Kp 
takes into account the dependence on the rigidity P of the GCR particle 
usually expressed in GV. In the present model Kp « P (e.g., see Ref.^°). 
Furthermore pk = 0.05 and, as described in Ref. 1, 



3B ' 



, , f 10 , m the polar regions, , , 

1(6') — i , . ^ . , . (2) 

1 , m the equatorial region. 

After the transformations from 3D field-aligned into 2D heliospherical 
coordinates^^, the symmetric components of the diffusion tensor contains 
both diagonal {Krr and Kee) and off-diagonal terms {Kr$ and Kgr), result- 
ing by a proper combination of K±0, K±r and K^^^. 



2. The Diffusion Paretmeter 

Kq accounts for the dependence on the solar activity. We estimated Kq by 
using the modulation strength (pg, in the framework of the Force Field (FF) 

approximation^^. (t>s was evaluated starting from Neutron Monitor (NM) 
counting rates in Ref. 13. We moreover used a practical correlation of Kq 
with the level of solar activity in the different solar phases^. We used, as 
solar activity monitor, the Smoothed Sunspot Number (SSN). 

This method is sensitive to the modulation of the GCR flux integrated 
over the full heliosphere, from the outer boundary to the Earth position, 
down to a lower limit in rigidity of ~ (2-3) GV. This limit is fixed by the 
sensitivity of the NM network, due to the geomagnetic rigidity cut-off and 
to the atmospheric yield function. The outer boundary of the heliosphere 
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is located at tho position of the Termination Shock. Beyond this hmit the 
model of heliosphere we are using is not more valid. Moreover, the additional 
modulation occurring in the heliosheat only affects particles with rigidity 
well bcilow 1 GV"'^'''"'^^. The method is sensitive also to the LIS used for the 
estimation of the modulation strength, but several LIS spectra do not differ 
each other above this rigidity limit. Finally the diffusion parameter depends 
on the outer boundary position, as follows from the FF approximation: 



l,^ Vs^{t) {Rts - -Rearth) r„s 
""'^'^ = • 

In Ref. 1 the boundary of heliosphere was placed at 100 AU. The solar 
cavity was split in 15 spherical regions to take into account the time spent 

by SW to travel outward. In each region of the interplanetary space, the 
parameters (i.e., SW speed, SSN, B^, tilt angle) are related to the solar 
activity at the time of the injection of the solar wind diffusing in that 
region-*^^. In this way modulated intensities of protons, down to ~ 400 MeV, 
were simulated and successfully compared with experimental data covering 
roughly one solar cycle. We did not find significant differences changing the 
position of the outer boundary of the heliosphere^. 



3. Heliospherical Size and Diffusion Parameter 

In the past years the position of the Termination Shock was estimated 
through the observations of Voyager 1 and Voyager 2 spacecrafts (see 
Refs. 17,18 and Table 1). In addition several authors (see Refs. 19,20) sug- 
gest that the size of the heliospere should change with the solar activity, 
following a quasi-periodic feature, roughly anti-correlated with the SSN. 



Table 1. Voyager crossings of Termination Shock. 





Rts (AU) 


solar latitude (deg) 


Voyager 1 


94.0 


+ 34.3 


Voyager 2 


83.7 


- 27.5 



Following these results we evaluated the effect of stretched and compressed 
heliospheres on the Cosmic Ray intensities at the Earth introducing a de- 
pendence of the diffusion parameter on the heliospherical size. We defined a 
new diffusion parameter Kq, introducing the parameter r(RTSi P) sensitive 
to the position of the Termination Shock: 
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K^{Rts) = r{RTS,P) i^o(100 AU) 



(4) 



r{RTs,P) = l + fiP) 



i?Ts (AU) - 100 
99 



(5) 



r{RTS,P) allows to modify the value of the diffusion parameter adapting it 
to a different volume of the heliosphere, determined by Rts- r{RTs, P) is 
fully effective below a rigidity limit Pi . We also defined a transition function 



For rigidity higher than P2, the dependence on Rts can be neglected. Here 
the diffusion parameter is still defined for an heliospherical dimension of 
100 AU. The dependence on the heliospherical radius Rts is then effective 
at rigidity lower than P2. Using the novel diffusion parameter Kq{Rts,P) 

wc simulated the modulated spectra, for different values of RtSj Pi and 
P2, extending the modulated spectra down to a lower rigidity. 

4. Results 

We compare our simulated spectra with proton data extended down to a ki- 
netic energy of 200 MeV. Here wc present results obtained using the follow- 
ing rigidity parameters: P2 = Pi = 1-0 GV. We used the Local Interstellar 
Spectrum (LIS) from Ref. 21 and compared it with the LIS form GAL- 
PROP^^. In Fig. 1 the results compared with AMS-01 data^^ are shown, 
assuming Rts = 120 AU, as discussed later on. 

We estimated the best value of Rts, looking at the RMS differences 
('7rms) with experimental data: 




for P > P2 , 
for Pi < P < P2 
for P < Pi . 



(6) 




with 



Vi = 




(8) 
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Fig. 1. Modulated proton spectra and eomparison with AMS-01 data'^'^. LIS are taken 
form Ref. 21 and GALPROP^'^. The vertical dotted line represents the lower limit of the 
sensitivity of the NM. Above this limit the two LIS are not significantly different. 



where Tj is the average energy of the i-th energy bin of the differential 
intensity distribution and an,i are the error bars including the experimental 
and Monte Carlo uncertainties. For each experimental spectrum we got the 
best values of Rts shown in Tabic 2 together with the minimum value of 
'7RMS- Data from BESS flights are given in Ref. 24, data from AMS-01 are 
given in Ref. 23. 

In Fig. 2 modulated spectra, obtained using values of Rts reported in 
Table 2, are shown in comparison with BESS experimental data. Modulated 
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spectrum compared with AMS-01 data has been shown in Fig. 1. We did 
not use data measured by Pamela^^ because pubhshed spectra start from 
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Table 2. Best values of Rts, its minimum and maximum values, and 
RMS differences between simulations and experimental data. 





(AU) 


R™^" (AU) 


(AU) 


'7RMS (%) 


BESS-1997 


115 


100 


130 


7.05 


AMS-1998 


120 


110 


135 


4.86 


BESS-1999 


120 


110 


130 


3.35 


BESS-2000 


140 


125 


150 


10.00 


BESS-2002 


105 


95 


115 


11.78 



400 MeV, while our analysis is more sensitive below this limit. In Tabic 2 wc 
report the interval of values of Rts where tjrms does not change by more 
than ~ (2 — 3) % from its minimum value, reported in the last column. This 
variation roughly represents the uncertainty of the computation itself, and 
it is determined comparing simulations and data at energies larger than 
(10 — 20) GeV, i.e. above the region of solar modulation. Results are shown 
in Fig. 3 in comparison with models^*^ and Voyager measurements^'"'^^. 

150 
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5' 110 
ct 100 
90 
80 
70 

1992 1994 1996 1998 2000 2002 2004 2006 2008 

Year 

Fig. 3. -Rts best value for the several experiments (crosses) in comparison with Voyager 
data^^'^* (dots) and models^" (dashed line). The shadow represents the region between 
the minimum and maximum value, as reported in Table 2. 
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As shown in Fig. 1 and Fig. 2, modulated spectra succeed to fit ob- 
serving data, in particular during periods of low solar activity. For more 
accurate results we need more accurate experimental data. Current error 
bars are of the order of 5% or even largcir. Moreover systematic deviations 
are present looking at different data sets at energy above (20 30) GV, where 
spectra are not affected by solar modulation. In addition current observa- 
tion data, taken on board of stratospheric balloons or space orbiters, may 
be contaminated at low energy by secondaries produced inside the Earth 
magnetosphere. A LIS spectrum with a slightly different shape could be also 
preferred to fit better low energy data. Finally a refinement of our model 
could be requested, starting with a slightly different values of P2 and Pi, in 
order to smooth the ripple present in some spectrum. In the future a model 
with an aspheric Heliosphere can be also developed. 

5. Conclusions 

We presented the HelMod 2D Monte Carlo code for the study of Cosmic 

Rays propagation in the inner heliosphere. Both heliospherical shape and 
size are supposed to be relevant for the modulation process. We introduced 
a dependence of the diffusion parameter on the heliospherical size, which 
accounts for the variation with time and solar activity. We compare mod- 
ulated spectra with experimental data covering the solar cycle 23. Then 
we found, for our 2D model, the best value of the heliospherical radius, 
which changes with time. Most of the solar modulation occurs in the inner 
heliosphere and differences in the heliospherical radius are effective only 
at energy below a few hundred MeV. Our results are not in contradiction 
with Voyager observations and models of TS distance as a function of solar 
activity. We found that LIS form Ref. 21 fits better observation data at low 
energy. 
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